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Abstract. Observations of the core of the massive cluster CI 0024+1654, at a redshift z ~ 0.39, were obtained with the Infrared 
Space Observatory using ISOCAM at 6.7 fjm (hereafter 7 yum) and 14.3 //m (hereafter 15 pm). Thirty five sources were detected 
at 15 yjm and thirteen of them are spectroscopically identified with cluster galaxies. The remaining sources consist of four stars, 
one quasar, one foreground galaxy, three background galaxies and thirteen sources with unknown redshift. The sources with 
unknown redshift are all likely to be background sources that are gravitationally lensed by the cluster. 

The spectral energy distributions (SEDs) of twelve cluster galaxies were fit from a selection of 20 models using the program 
GRASIL. The ISOCAM sources have best-fit SEDs typical of spiral or starburst models observed 1 Gyr after the main starburst 
event. The star formation rates were obtained for cluster members. The median infrared luminosity of the twelve cluster galaxies 
is ~ 1.0 x 10" L G , with 10 having infrared luminosity above 9 x 10 10 L Q , and so lying near or above the 1 x 10" L G threshold 
for identification as a luminous infrared galaxy (LIRG). The [O n] star formation rates obtained for 3 cluster galaxies are one to 
two orders of magnitude lower than the infrared values, implying that most of the star formation is missed in the optical because 
it is enshrouded by dust in the starburst galaxy. 

The cluster galaxies in general are spatially more concentrated than those detected at 15 fim. However the velocity distribu- 
tions of the two categories are comparable. The colour-magnitude diagramme is given for the galaxies within the ISOCAM 
map. Only 20% of the galaxies that are significantly bluer than the cluster main sequence were detected at 15 fjm, to the limit- 
ing sensitivity recorded. The counterparts of about half of the 15 /jm cluster sources are blue, luminous, star-forming systems 
and the type of galaxy that is usually associated with the Butcher-Oemler effect. HST images of these galaxies reveal a dis- 
turbed morphology with a tendency for an absence of nearby companions. Surprisingly the counterparts of the remaining 15 yum 
cluster galaxies lie on the main sequence of the colour-magnitude diagramme. However in HST images they all have nearby 
companions and appear to be involved in interactions and mergers. Dust obscuration may be a major cause of the 15 fjm sources 
appearing on the cluster main sequence. The majority of the ISOCAM sources in the Butcher-Oemler region of the colour- 
magnitude diagram are best fit by spiral-type SEDs whereas post-starburst models are preferred on the main sequence, with the 
starburst event probably triggered by interaction with one or more galaxies. 

Finally, the mid-infrared results on CI 0024+1654 are compared with four other clusters observed with ISOCAM. Scaling 
the LIRG count in CI 0024+1654 to the clusters Abell 370, Abell 1689, Abell 2218 and Abell 2390 with reference to their virial 
radii, masses, distances, and the sky area scanned in each case, we compared the number of LIRGs observed in each cluster. 
The number in Abell 370 is smaller than expected by about an order of magnitude, even though the two clusters are very similar 
in mass, redshift and optical richness. The number of LIRGs detected in each of Abell 1689, Abell 2218 and Abell 2390 is 0, 
whereas 3 were expected from the comparison with CI 0024+1654. A comparison of the mid-infrared sources in Abell 1689 
and Abell 2218 shows that the sources in Abell 1689 are more luminous and follow the same trend identified in the comparison 
between CI 0024+1654 and Abell 370. These trends seem to be related to the dynamical status and history of the clusters. 

Key words. Galaxies: clusters: general - Galaxies: clusters: individual (CI 0024+1654) - Infrared: galaxies 



1. Introduction 

Clusters of galaxies contain thousands of members within a re- 
gion a few Mpc in diameter, and are the largest known grav- 
itationally bound systems of galaxies, having masses up to 
10 15 M for the richest systems. In hierarchical models clus- 
ters of galaxies grow by accreting less massive groups falling 
along filaments at a rate governed by the initial density fluc- 
tuation spectrum, the cosmological parameters and the nature 
and amount of dark matter. In the cluster environment, newly 
added galaxies are transformed from blue, active star forming 
systems, to red, passive ellipticals, undergoing a morphological 
evolution stronger than that of field galaxies at a similar redshift 
(Gavazzi & Jaffe ll987l Byrd & Valtonen [T990l Abraham et al. 
1996b). The cluster galaxy population is also characterized by 
a lower star formation rate (SFR) than field galaxies of simi- 
lar physical size and redshift (Couch et al. 120011 Lewis et al. 

nom 

Butcher and Oemler ( 1978) showed that clusters of galax- 
ies generally have a fraction fs of blue galaxies 1 that increases 
with cluster redshift, ranging from a value near at z — 0, to 
20% at z = 0.4 and to 80% at z = 0.9, suggesting a strong 
evolution in clusters (Rakos & Schombert 119951 . The galax- 
ies responsible for the Butcher-Oemler effect (hereafter BO 
effect) are generally luminous, spirals, and emission-line sys- 
tems, with disturbed morphologies. The high resolution imag- 
ing achieved by the Hubble Space Telescope (HST) has greatly 
improved morphological studies of the cluster galaxy popu- 
lation over a wide range in redshift. The mixture of Hubble 
types in distant clusters is significantly different from that seen 
in nearby systems. The population of star-forming and post- 
starburst galaxies are disk dominated systems, some of which 
are involved in interactions and mergers. (Abraham et al 1996a 
Stanford et al. [79981 Couch et al. [T9981 Van Dokkum et al. 
rH?9lfl Morris et al. H^gl Poggianti et al. [T9^9l Best l2U00l 

Many mechanisms have been proposed to explain the com- 
plicated processes that occur in clusters, including ram pres- 
sure stripping of gas (Gunn & Gott I1972> . galaxy harass- 
ment (Moore et al. 1 19961 Moss & Whittle [1997), galaxy infall 
(Ellingson et al. 2001 ), cluster tidal forces (Byrd & Valtonen 
119901 Fujita[T998 1 and interactions with other cluster galaxies 
(Icke 1985 Moss & Whittle lT997l . The main processes respon- 
sible for the morphological and spectral evolution of cluster 
galaxies have yet to be determined. Ram pressure and tidal ef- 
fects can quench the star formation activity gradually because 
they operate over a period longer than 1 Gyr (e.g. Ghigna et al. 
119981 Ramirez & de Souza 1998). Galaxy-galaxy interactions, 
galaxy harassment and cluster mergers can enhance it and pro- 
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1 Blue galaxies are defined as brighter than M v = -19.26 (with 
Ho = 70 km s~' Mpc -1 ) with rest-frame B-V colours at least 0.2 mag- 
nitudes bluer than those of the E/S0 galaxy sequence at the same ab- 
solute magnitude (Butcher & OemlerfTWJ Oemler et al. 119971 . 



duce changes in galaxy properties over timescales of ~ 100 
Myr (e.g. Lavery & Henry 119861 Moore et al. I1996i . Recent 
changes in the properties of the galaxies may be detectable in 
the mid-infrared if associated with a burst of star formation. 

In the context of our ongoing exploitation of mid-infrared 
cluster data obtained with ESA's Infrared Space Observatory 
(ISO, Kessler et al. 1996), we have analysed ISO observations 
of Abell 370 , Abel l 2218 & Abel l 2390 (Metcalfe et al. l2003l 
Altieri et al. fT999l Biviano et al. l2004l Abell 2219 (Coia et 
al. I2004> . and CI 0024+1654 (this paper). In common with 
other surveys we have found that mid-infrared observations of 
field galaxies from deep surveys reveal a population of starburst 
galaxies that evolve significantly with redshift (Aussel et al. 
IT9991 Oliver et al. l2000l Serjeant et al. l2000l Lari et al. l200Tl 
Gruppioni et al. l200l Elbaz & Cesarskv 120031 Metcalfe et al. 
2003 1 Sato et al. 2003 ). This class of sources are Luminous and 
Ultraluminous infrared galaxies (LIRGs and ULIRGs, Sanders 
& Mirabel [T9961 Genzel & Cesarsky 120001 have SFRs of 
~ 100 M G yr 1 (Oliver et al. l2000l Mann et al. l2002T ) and seem 
to be almost always the result of galaxy-galaxy interactions at 
least in the local Universe (Veilleux et al. l2002l ). The ISOCAM 
sources account for most of the contribution of the mid-infrared 
to the Cosmic Infrared Background (CIRB, Altieri et al l 19991 
Franceschini et al. l200Tl Metcalfe et al. l20TTTll2003l Elbaz et al. 
2002). Studies of the global SFR show a decline by a factor of 
3-10 since the peak of star formation at z - 1-2 (Madau et al. 
fT996l Steidel et al. lT^l Elbaz & Cesarskv I2UU31 . The down- 
turn in the global SFR and population of LIRGs and ULIRGs 
may be caused by galaxies running out of gas available for 
star formation and the buildup of large scale structure in the 
Universe that changed the environment of galaxies. 

A study of the impact of the environment on galaxies in 
clusters could help in understanding the global SFR. Mid- 
infrared observations have been published for local clusters 
(Boselli et al. fT997lfT998l C onturs i et al. l200Tl an d dista nt clus- 
ters of galaxies (Pierre et al. 1996 Lemonon et al. 1998 Altieri 
et al. fT999l Fadda et al. l2000l Metcalfe et al. l2003l Coia et al. 
2004}. In Abell 2390, Abell 370 and Abell 2218 the greater 
part of the 15 fim sources for which spectroscopic redshifts are 
available are found to be background sources (Metcalfe et al. 
120031 However Fadda et al. ( 2000 ) and Due et al. ( 12002b found 
a higher proportion of 15 fim cluster sources in the cluster Abell 
1689 at z = 0.18 and Due et al. (2004) discovered many LIRGs 
in the cluster J1888.16CL at z = 0.56. 

In this work we focus on the mid-infrared properties of the 
galaxy cluster CI 0024+1654. The paper is organized as fol- 
lows: Sect. |21 contains a description of the cluster. Section |5] 
describes the infrared observations and outlines the data re- 
duction, source extraction and calibration processes. Section^] 
presents the results, the model spectral energy distributions 
(SEDs) and star formation rates for cluster galaxies. Section|5] 
describes the spatial, redshift and colour properties of cluster 
galaxies and contains a description of Hubble Space Telescope 
images of some galaxies detected by ISOCAM. Section [6] 
makes a comparison between CI 0024+1654 and other clus- 
ters studied, including Abell 1689, Abell 370, Abell 2390 and 
Abell 2218. The conclusions are in Sect.0 The Appendix con- 
tains additional comments on some of the ISOCAM sources. 
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Fig. 1. The 15 yum contour map (red) of CI 0024+1654 overlaid on a Very Large Telescope image taken in the V band with the 
FORS2 instrument (ESO program identification: 65.O-0489(A)). Numbers 1 to 30 refer to the 15/mi sources in the primary 
list (Table |2j and sources Is to 5s label sources from the supplementary list (Table [3j. Each set of sources is labelled in order 
of increasing Right Ascension. Sources ISO_C10024_29 and ISO_C10024_30 are outside the boundary of the optical map, and 
have a star and a faint galaxy, respectively, as optical counterparts. Blue circles denote 15/mi sources that are spectroscopically 
confirmed cluster galaxies. Greek letters (a + 6) identify four gravitationally lensed images of the background galaxy associated 
with the spectacular giant arcs in the cluster field. North is up and East is to the left. The centre of the ISO map, indicated by a 
cross, is at R.A. 00 26 37.5 and DEC. 17 09 43.4 (J2000). 



We adopt H = 70 km s" 1 Mpc 1 , Q A = 0.7 and fi m = 0.3. 
With this cosmology, the luminosity distance to the cluster is 
Dl = 2140 Mpc and 1" corresponds to 5.3 kpc at the cluster 
redshift. The age of the Universe at the cluster redshift of 0.39 
is 9.3 Gyr. 

2. The cluster 

Cl 0024+1654 is a rich cluster of galaxies at redshift z ~ 0.395 
(Humason & Sandage [T957l Gunn & Oke [T975l Smail et al. 
119931 . It has a spectacular system of gravitationally lensed arcs 
(Fig. [0 that were first observed by Koo (1988 1 and subse- 
quently studied by Mellier et al. ( H991L Kassiola et al. ©94 1, 
Wallington et al. (1995) and Smail et al. (19971. The main arc 
is split into segments (Colley et al. 119961 Tyson et al. Q998 1 
and is the lensed image of a background blue galaxy at redshift 



Z = 1.675 (Broadhurst et al. 2000). Unlike other clusters of 
galaxies, such as Abell 2218 or Abell 2390 that have arc-like 
features, Cl 0024+1654 does not have a central dominant cD 
galaxy. 

Cl 0024+1654 was one of the two clusters examined by 
Butcher & Oemler J1978b in their first published work on 
galaxy colours. The fraction of blue galaxies in Cl 0024+1654 
is fe = 0.16 which is much larger than the values of fs = 0.04 
and f B = 0.03 for the Virgo (z = 0.003) and Coma (z = 0.02) 
clusters respectively (Butcher & Oemler 1984 Dressier et al. 
IT9851 Schneider et al. fT986l . 

The determination of the spectral types of galaxies has been 
made for the cluster and the field (Czoske et al. 2002 Balogh 
et al. l2002l. Usin g a wide field HST survey of Cl 0024+1654, 
Treu et al. ( 2003 1 found that the fraction of early-type galaxies 
(E+S0) is highest (~ 73%) in the cluster core, declines rapidly 
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Fig. 2. The fractional distribution of velocities of cluster galax- 
ies in CI 0024+1654. The counterparts of the cluster members 
detected at 15 /mi are displayed in bold. All velocities are in the 
cluster rest-frame. 



in redshift of the two main galaxy concentrations, is about 
3000 kms -1 . The collision scenario is not unusual because 
30-40% of galaxy clusters have substructures that are detected 
in optical, X-ray and radio observations (e.g. Dressier & 
Shectman lT988l Schuecker et al. 1200 II Mercurio et al. l2003> . 
These substructures indicate that the hosting clusters are not 
fully dynamically relaxed and might be currently undergoing, 
or have recently undergone, mergers. The important results 
obtained by Czoske et al. (2001 , 2002) may give a new insight 
into the BO effect in medium and high redshift clusters. The 
impact of the two clusters could have triggered star formation 
in cluster galaxies, especially those on the leading edge of 
the smaller cluster, thus generating the large number of blue 
galaxies observed in CI 0024+1654. The differences between 
the results from the methods used to determine the mass of 
CI 0024+1654 are naturally explained by the cluster collision 
(Czoske et al. lSUOlb . 



to about 50% at ~ 1 Mpc, and reaches the background value of 
~ 43% at the periphery of the cluster (at a radius of ~ 5 Mpc). 

The mass profile of the cluster has been inferred from grav- 
itational lensing analyses, kinematical analyses of the redshifts 
of the cluster galaxies, and X-ray observations. The low value 
of the X-ray luminosity implies a mass 2 to 3 times smaller than 
that predicted by lensing models (Soucail et al. 2000 Ota et al. 
2004 1. There are also differences among the lensing models. 
Broadhurst et al. {2000 1 favour a cuspy NFW (Navarro, Frenk 
& White 1997 1 profile for the mass distribution of the clus- 
ter. Tyson et al. ©98 1 argue that this result requires an aver- 
age cluster velocity distribution much higher than the measured 
value of cr v ~ 1 150km s _1 (Dressier et al. 1999 > and favour a 
uniform-density core for the halo mass profile. The weak lens- 
ing analysis made by Kneib et al. ( 2003 ), based on an extensive 
HST survey of CI 0024+1654 (Treu et al. l2003l . indicates that 
the region within 5 Mpc from the centre of the cluster is well fit 
by a steep NFW-like profile, and the isothermal fit is strongly 
rejected. Ota et al. ( 2004 1 analysed the X-ray emission from the 
cluster and found that an isothermal model with a mean tem- 
perature of 4.4 keV is a good fit to the temperature distribution. 
They also estimated the mass of the cluster on the assumption 
that the intracluster medium is in hydrostatic equilibrium, and 
confirmed that the mass inferred from X-ray observations is 
smaller than the mass predicted by lensing models by a factor 
of about 3. The complex structure of the central region of the 
cluster might explain the discrepancy (Zhang et al. 2004). 

During an extensive spectroscopic survey of CI 
0024+1654, Czoske et al. J200H detected a significant 
number of galaxies at a redshift slightly lower than that of 
the cluster (Fig. |3. This concentration (component B) is 
interpreted as a small, less massive cluster superimposed on 
the main cluster (component A) and lying at a mean redshift of 
Z = 0.381. Numerical simulations indicate that the two clusters 
in CI 0024+1654 are involved in a high-speed collision along 
the line-of-sight to the cluster (Czoske et al. 2002 (. The relative 
velocity of the two components, derived from the difference 



3. Observations, data reduction and source 
detection 

3.1. Observations 

The core of the cluster CI 0024+1654 was observed at 7 fim and 
15 /im using the LW2 and LW3 filters of the camera ISOCAM 
(Cesarsky et al. [T996l ) on board ISO. The ISOCAM long wave- 
length (LW) detector consisted of a 32 x 32 SiGa array. The 
observations were made on June 9 and June 15, 1997, in raster 
mode, with an on-chip integration time of 5.04 seconds in the 
3" per pixel field of view, and cover an area of approximately 
38 square arcminutes. A discussion of the relevant observing 
strategy for ISOCAM can be found in Metcalfe et al. (2003) 
(with the exception that the CI 0024+1654 rasters were not 
"micro-scanned" with sub-pixel finesse. I.e. the raster step sizes 
were multiples of the array pixel size, whereas Metcalfe et al. 
describe rasters with step sizes involving fractions of a pixel's 
dimensions). For each pointing of the raster, 14 readouts were 
performed at 7 fim and 10 at 15 fim. Fifty readouts were taken at 
the beginning of each raster to allow for detector stabilization. 
The parameters used for the observations are given in Tabled 
which also lists the 5<x sensitivity reached after data processing. 
The diameter of the point spread function (PSF) central max- 
imum at the first Airy minimum is 0.84 x A{fim) arcseconds. 
The FWHM is about half that amount and Okumura ([1998) 
obtained values of 3.3" at 7yum and 5" at 15jt/m for the PSF 
FWHM in the 3" per pixel field-of-view. The data were reduced 
and, to take advantage of a slight deliberate offset between the 
two 15 fim maps, were rebinned so that the final map has a pixel 
size of 1", with potential for slightly improved spatial resolu- 
tion. The maximum depth is reached toward the centre of the 
rasters, where the dwell time per position on the sky is great- 
est. As recorded in Tabled two rasters were made at 15/im 
and a single raster at 7 fxm that had a larger step size and less 
observation time (by a factor of 7). 
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Table 1. Observational parameters used. The observations were made with the ISOCAM LW2 (5 - 8.5 /mi) and LW3 (12 - 18 //m) 
filters at their respective reference wavelengths of 6.7 /mi and 14.3 /urn. On-chip integration time was always 5.04 seconds and the 
3" per-pixel-field-of-view was used. M and N are the number of steps along each dimension of the raster, while dm and dn are the 
increments for a raster step. The sensitivity in //Jy is given at the 5cr level. The table also includes the total area covered and the 
number of readouts per raster step. The rasters were repeated k times. Tot. t is the total time dedicated to each filter observation. 



Filter 


'iref 


n Steps 


dm 


dn 


Reads 


area 


Done k 


Sensitivity 


Tot. t 




dim) 


M N 


(") 


(") 


per step 


C 2 ) 


times 


0*Jy) 


(sec) 


LW2 


6.7 


6 6 


45 


45 


14 


28.6 


1 


400 


3138 


LW3 


14.3 


14 14 


21 


21 


10 


37.8 


2 


140 


22615 



3.2. Data reduction 

The data were reduced using the ISOCAM Interactive Analysis 
System or CIA (Delaney & Ott l2002l Ott et al. Il997l in con- 
junction with dedicated routines, following the method 2 de- 
scribed in Metcalfe et al. (|2p03l. The two 15 /mi raster maps 
were normalized to their respective redundancy maps 3 and 
merged into a single raster, thus increasing the sensitivity of 
the map to faint sources. 

The 15 /mi sources were found by visually inspecting the 
merged and the two individual maps. Detections common to 
the three maps are given in the primary source list (Table 0. 
The selection criterion adopted is very conservative because it 
requires the source to be detected in both individual rasters, 
resulting in the rejection of faint sources that can be detected 
only in the merged map. Therefore we include a supplemen- 
tary list containing sources found in the merged map and in at 
least one of the individual rasters (Table QJ. The fluxes of the 
sources were computed by aperture photometry using a circular 
aperture of 9" diameter, centred on the infrared source. 

Since there was only one raster at 7/im, an independent 
comparison of source detections was not possible and real 
sources were considered to be only those having a 15 /mi coun- 
terpart. 

The optical data were obtained on the ESO/VLT dur- 
ing October 2001 with the VLT/FORS instrument in service 
mode 4 . The data have been processed at the TERAPIX data 
center 5 . Pre-calibrations, astrometric and photometric calibra- 
tions as well as image stacking and catalog production were 
done using standard CCD image processing algorithms and 

2 An interesting alternative approach, optimised for extended- 
source detection, but which gives very good overall results and can 
be used to cross-check the validity of the faintest detections, is the 
method of Miville-Deschenes et al. <2000t . We employed products of 
this alternative reduction process (SLICE) as a further cross-check of 
the reality of the faintest sources in our lists, reasoning that both meth- 
ods would tend to preserve real sources, while any residual glitches af- 
fecting the results of our primary analysis would have some probabil- 
ity of being filtered out in the independent SLICE analysis. To arrive 
at our final list we rejected some very faint candidate sources as be- 
ing unreliable when they were poorly reproduced by this independent 
analysis route. We consider this approach to be very conservative. 

3 The "redundancy" of a point in a raster map refers to the number 
of raster steps for which that point on the sky has been seen by some 
detector pixel. 

4 Program ID: 65.0-489A FORS; PI: Fort 
3 http://terapix.iap.fr 



tools available at the TERAPIX center (see McCracken et al. 
2003 for details). 

3.3. Monte Carlo simulations and calibration 

Monte Carlo simulations were performed to calibrate the com- 
plex data reduction process by characterizing how it affects 
model point sources with known properties inserted into the 
raw data, and to characterize the precision with which source 
signals could be recovered from the data. The procedure 
adopted is fully described in Metcalfe et al. ( 2003 1, the only 
exception being that all photometry performed for the present 
work employed the XPHOT routine in CIA, whereas Metcalfe 
et al. used SExtractor (Bertin & Arnouts 119961 . The fluxes 
of the inserted fake sources were measured in exactly the 
same way as the fluxes of the real sources. The simulations 
were performed independently for all rasters. The fluxes of 
the inserted fake sources were measured by manual photom- 
etry adopting the same parameters, i.e. same aperture, used for 
the real sources. The simulations establish the relationship be- 
tween source signal detected in the photometric aperture, and 
total source signal actually collected in the detector. Signal in 
the photometric aperture can therefore be scaled to total signal, 
but at this point the units are detector units, or ADUgs. The re- 
lationship between source signal deposited in the detector and 
actual source flux-density in mJy is determined by two further 
scaling factors, (a) A filter specific ISOCAM calibration factor 
relating stabilised source signal, in detector units, to mJy, i.e. 1 
ADU per gain per second (ADUgs) = 0.43 mJy for LW2 and 1 
ADUgs ee 0.51 mJy for LW3 (Delaney & Ott l2002l : and (b) a 
correction for detector responsive transient effects which cause 
the signal from faint sources to fall below that expected on the 
basis of bright reference source measurements. The derivation 
of that scaling factor, in the context of the reduction algorithm 
applied here, is described in Metcalfe et al. ( 2003 ) and refer- 
ences therein. 

4. Results 

The merged 15 /mi map (Fig.^ is overlaid on a V-band image 
of the cluster taken with the Very Large Telescope (VLT). 

The list of sources detected in the two individual 15//m 
maps is given in Table [2] and a supplementary list is given in 
Table|5](see Sect.[3J. The name of the ISOCAM source is com- 
posed of the satellite acronym (ISO), the partial name of the 
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Table 2. The primary list of sources detected at 15 fim. The columns are (from left to right): source identification number in 
order of increasing Right Ascension; source flux and precision before calibration (in ADUgs); source flux and precision after 
calibration (in mJy); Right Ascension and Declination of the infrared source (J2000); redshift and name of the optical counterpart, 
if known. The redshift of source ISO_C10024_02 was taken from Schmidt et al. J1986I . All other redshifts are from Czoske et al. 
J200U or provided by T. Treu (private communication). The names of the stars are from HEASARC. 



Source ID 


Signal 


Precision 


Flux 


Precision 


R.A. 


DEC. 


z 


Name of optical 


ISCLC10024 


(ADU) 


(± ADU) 


(mJy) 


(± mJy) 


(J2000) 


(J2000) 


if known 


counterpart if known 


01 


0.634 


0.090 


0.786 


0.120 


00 26 24.3 


+ 17 08 18.0 


Star 


TYC 1180-82-1 


02 


0.472 


0.060 


0.577 


0.090 


00 26 26.3 


+ 17 09 39.0 


0.959 


PC 0023+1653 


03 


0.256 


0.060 


0.298 


0.080 


00 26 27.0 


+ 17 09 11.7 


Star 


N323 1320329 


04 


0.298 


0.060 


0.352 


0.080 


00 26 30.8 


+ 17 09 44.8 


_ 




05 


0.184 


0.060 


0.205 


0.080 


00 26 30.7 


+ 17 09 03.5 


0.3935 




06 


0.230 


0.060 


0.264 


0.080 


00 26 30.9 


+ 17 10 29.4 


_ 




07 


0.217 


0.060 


0.247 


0.080 


00 26 31.0 


+ 17 09 37.4 


_ 




08 


0.610 


0.090 


0.755 


0.120 


00 26 31.0 


+ 17 10 15.6 


0.4005 




09 


0.286 


0.060 


0.336 


0.080 


00 26 31.5 


+ 17 1021.4 


0.2132 




10 


0.191 


0.060 


0.214 


0.080 


00 26 31.8 


+ 17 10 03.5 


0.4000 




11 


0.220 


0.060 


0.251 


0.080 


00 26 35.5 


+ 17 09 59.1 


0.5558 




12 


0.182 


0.060 


0.202 


0.080 


00 26 35.9 


+ 17 10 36.4 


0.3860 




13 


0.210 


0.060 


0.238 


0.080 


00 26 35.9 


+ 17 11 39.9 






14 


0.293 


0.060 


0.346 


0.080 


00 26 36.6 


+ 17 07 44.4 


Star 


EO903-0219757 


15 


0.170 


0.060 


0.186 


0.080 


00 26 36.6 


+ 17 10 25.4 


_ 




16 


0.260 


0.060 


0.302 


0.080 


00 26 37.7 


+ 17 11 33.4 






17 


0.253 


0.060 


0.293 


0.080 


00 26 38.9 


+ 17 08 03.9 


_ 




18 


0.634 


0.090 


0.786 


0.120 


00 26 40.4 


+ 17 09 28.7 


0.394 




19 


0.187 


0.060 


0.208 


0.080 


00 26 40.6 


+ 17 11 44.9 


. 




20 


0.262 


0.060 


0.305 


0.080 


00 26 40.7 


+ 17 09 54.0 


0.7125 




21 


0.191 


0.060 


0.213 


0.080 


00 26 41.3 


+ 17 11 01.4 


0.3924 




22 


0.334 


0.060 


0.398 


0.080 


00 26 41.8 


+ 17 09 52.9 


0.3935 




23 


0.282 


0.060 


0.331 


0.080 


00 26 42.7 


+ 17 09 41.6 






24 


0.189 


0.060 


0.211 


0.080 


00 26 42.8 


+ 17 08 48.5 


0.3954 




25 


0.200 


0.060 


0.225 


0.080 


00 26 42.9 


+ 17 08 31.7 


0.9174 




26 


0.324 


0.060 


0.385 


0.080 


00 26 43.1 


+ 17 08 25.1 


0.3961 




27 


0.373 


0.060 


0.449 


0.080 


00 26 43.2 


+ 17 09 11.7 


0.3932 




28 


0.120 


0.060 


0.122 


0.080 


00 26 45.5 


+ 17 08 30.8 






29 


0.946 


0.110 


1.191 


0.140 


00 26 47.9 


+ 17 08 12.0 


Star 


EO903-0219375 


30 


0.489 


0.070 


0.598 


0.090 


00 26 48.0 


+ 17 09 15.7 







Table 3. The supplementary list of 15 fim sources. The meaning of the columns is the same as in Table[2] 



Source ID 


Signal 


Precision 


Flux 


Precision 


R.A. 


DEC. 


z 


ISOX10024 


(ADU) 


(± ADU) 


(mJy) 


(± mJy) 


(J2000) 


(J2000) 


if known 


Is 


0.150 


0.060 


0.160 


0.080 


00 26 31.6 


+ 17 09 17.1 


0.3998 


2s 


0.190 


0.060 


0.212 


0.080 


00 26 32.4 


+ 17 09 06.5 




3s 


0.124 


0.060 


0.127 


0.080 


00 26 32.6 


+ 17 10 26.1 


0.3965 


4s 


0.173 


0.060 


0.190 


0.080 


00 26 37.4 


+ 17 09 29.1 


0.3900 


5s 


0.182 


0.060 


0.202 


0.080 


00 26 42.9 


+ 17 10 45.6 





Table 4. The list of 7 /im sources with source flux exceeding 5<x of the local noise, and having 15 jim counterparts. From left to 
right: 15 /urn source identification number as listed in Table |2j signal and precision of flux in ADU; flux and precision in mJy; 
Right Ascension, Declination; nature of the source; and finally the [7 yum]/[15 jum] colour ratio. 



Source ID 


Signal 


Precision 


Flux 


Precision 


R.A. 


DEC. 


Notes 


[7^m]/[15^m] 




(ADU) 


(± ADU) 


(mJy) 


(± mJy) 


(J2000) 


(J2000) 






03 


1.461 


0.219 


1.685 


0.253 


00 26 27.0 


+ 17 09 08.2 


Star 


5.7 


14 


1.276 


0.191 


1.468 


0.220 


00 26 36.5 


+ 17 07 41.0 


Star 


4.2 


29 


3.806 


0.571 


4.427 


0.664 


00 26 47.9 


+ 17 08 11.1 


Star 


3.7 
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Table 5. References to ISOCAM sources in existing archives. 
The first column lists the source number from Tables |2] and [3] 
The key to the references in Cols. 2 to 8 is Dressier & Gunn 
fT992l DG), Czoske et al. d2"0UTl CKS), Schneider et al. ( fT986l 
SDG), Smail et al. STW7\ SDC), McLean & Teplitz ( fT9%l 
MT), Butcher & Oemler JT9781 BO ), Soucail et al. J2UUD1 S) 
and Pickles & van der Kruit ( 119911 P). The table does not in- 
clude the four stars or the quasar PC 0023 + 1653. 



id 


DG 


CKS 


SDG 


SDC 


MT 


BO 


Others 


05 


198 


262 


223 






123 




06 


278 








64 






08 


264 


267 


113 




22 


34 


S01 


09 


257 


282 






6 






10 


237 


289 












12 




381 












17 


12 














18 


47 


444 




797 


4 






21 


119 


453 












22 


43 


459 


146 


834 








23 


23 






883 








24 




471 












27 


3 


474 




928 








Is 


195 


280 






50 






3 s 


246 


304 






55 






4s 


101 


412 


186 


573 


11 


87 


P61 



cluster (C10024), and an identification number assigned to each 
source. 

The 7 yum fluxes for the stars, and their [7 /im]/[15 /mi] flux 
ratios, are given in Table 0] The labels of the 7 jjm sources are 
taken to be the same as their 15 fim counterparts. The fourth star 
detected at 15 pm (ISO_C10024_01) is outside the boundary of 
the 7 yi/m raster. No other sources were detected above the 5 cr 
limitof ~400yuJy. 

The column named "Precision" in Tables |2] and [3] reflects 
the repeatability of the photometric results in the fake source 
simulations. For each source brightness the precision is the lcr 
scatter found in the recovered fluxes of fake sources of similar 
brightness inserted into the raw data. It should be noted that the 
ratio signal/precision is not, however, the signal-to-noise ratio; 
nor is it the significance of a source detection. The precision, as 
defined here, includes, for example, signal deviations caused by 
residual glitches around the map and by variations in residual 
background level across the map, which are of course factors 
that can influence the accuracy of source photometric measure- 
ments. However, the significance of an individual source detec- 
tion must be determined through an examination of the pixel- 
to-pixel signal variations in the map local to the source. All of 
the sources listed in Tables [2] and [3] have significance values 
greater than 5<x. The quoted precision is a measure of the rela- 
tive calibration accuracy within the sample recorded. The abso- 
lute flux calibration is, in addition, affected by factors such as 
imperfections in the knowledge of the detector response and re- 
sponsive transient correction, and differences between the spec- 
tral shapes of measured sources with respect to the canonical 
calibration source for ISOCAM, which has a stellar spectrum. 
Experience over many surveys suggests an absolute calibration 



accuracy of + 15% for ISOCAM faint-source photometry due 
to accumulated systematic effects. The measured colour ratios 
for the stars listed in Table |4] are consistent with this, and give 
some confidence in the absolute calibration at the ±15% level. 

The number of sources in the primary and supplementary 
lists are 30 and 5, respectively, yielding a total of 35 sources. 
A search was performed for counterparts of the ISOCAM 
sources at other wavelengths, using a search radius of 6" in 
HEASARC 6 and NASA/IPAC Extragalactic Database 7 . Of the 
30 sources in the primary list (Table [2J|, 19 are identified with 
catalogued sources. Of these nineteen, four are stars, one is a 
quasar (PC 0023+1653), ten are cluster members, one is a fore- 
ground galaxy and three are background galaxies. As for the 5 
supplementary sources (Tabled, three are identified with clus- 
ter galaxies. There is a total of 13 sources with unknown red- 
shift and all have optical counterparts (Fig.Q. 

The median flux-density of the unknown-redshift sources is 
about 250//Jy. The expected number of 15 fim sources lensed 
by the cluster down to that flux limit is about 20 ± 10, based on 
the log N - log S distribution of Metcalfe et al. (2003 ). So the 
expected number of background sources is sufficient to account 
for the observed number of unknown-redshift sources. 

4.1. Spectral Energy Distributions 

Spectral energy distributions were computed for the cluster 
galaxies using the 15/um fluxes in Tables |2 and [3] and a 
5<x upper limit of 400 fiJy for the 7 fim fluxes. The spec- 
tral range of the ISOCAM sources was extended by includ- 
ing measurements in the optical photometric bands and the 
near-infrared. The values were retrieved from the NASA/IPAC 
Extragalactic Database. The references to extensive observa- 
tions of the ISOCAM sources from a wide range of catalogues 
are listed in Table [5] 

The SEDs were modelled using the program GRASIL 
(Silva et al. I1998> . These models have already been used by 
Mann et al. ( I2002> to fit the SEDs of galaxies detected by 
ISOCAM in the Hubble Deep Field South. Given the lim- 
ited photometric accuracy and, for some galaxies, the limited 
amount of available photometric data, it is not possible to ob- 
tain an accurate model of the SED of each individual galaxy. 
Instead, observed SEDs are compared with models representa- 
tive of broad classes of spectral types. The number of models 
to be considered needs to be adequate to reflect the quality of 
the data and the questions to be answered (see, e.g., Rowan- 
Robinson 2003 (. The redshifts of the galaxies are known and 
hence it was possible to use more than the five models of Mann 
et al. (2002 1 without the risk of multiple solutions. 

In this work, 20 models were considered. They were taken 
from either the public GRASIL library, or built by running 
the publicly available GRASIL code, or kindly provided by 
L. Silva (private comm.). The 20 models reproduce the SEDs 
of several kinds of galaxies, from early-type, passively evolv- 

6 http://heasarc.gsfc.nasa.gov/db-perl/W3Browse/w3browse.pl 

7 The NASA/IPAC Extragalactic Database (NED) is operated by 
the Jet Propulsion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space Administration. 
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Fig. 3. SEDs for cluster members in CI 0024+1654 detected at 15 /im. The horizontal axis is the wavelength in the cluster rest 
frame and the vertical axis is the flux density (in normalized units). The data points are given by black dots and the model fit to 
the SED by a continuous line. 



ing ellipticals (labelled 'E', in the following), to spiral galaxies 
(labelled 'S'), and starburst galaxies (similar to the local ex- 
amples Arp220 and M82), as seen either at the epoch of the 
starburst event, or 1 Gyr after (labelled 'SB' or 'SB+1', respec- 
tively). At variance with Mann et al. (2002i we avoid consid- 
ering models older than the age of the Universe at the cluster 
redshift. Specifically, we consider models for galaxies with a 
formation redshift z = 1 (this is the redshift when stars start 
forming in the GRASIL models), corresponding to an age at 
the cluster redshift of ~ 3 Gyr, and models with a formation 
redshift z = 4, corresponding to an age at the cluster redshift of 
~ 8 Gyr. 

The best-fitting SED model was determined by x 1 mini- 
mization, leaving the normalization of the model free. In prin- 
ciple, a given SED model cannot be freely rescaled, since the 
parameters in the GRASIL code depend on the mass of the 
galaxy. However, a fine tuning of the values of the GRASIL 
parameters only makes sense if the photometric uncertainties 



of the observed galaxy are small enough, and the SED is very 
well covered by observations, which is not the case here. As a 
matter of fact, the model SEDs of two galaxies of very different 
mass, such as M82 and NGC6090, could not be distinguished 
after rescaling, in the wavelength range covered by our obser- 
vations, and with the accuracy of our photometric data points. 

The SEDs are plotted in Fig. [3] and summarized in Table 
[6] for all cluster galaxies detected with ISOCAM, except for 
ISO_C10024_12, for which there are not enough photometric 
data. All the fits are acceptable at the 1 % confidence level, and 
all, except two, at the 5% confidence level. S and SB+1 models 
provide the best-fit to most SEDs, but the fits are, in general, not 
unique. Among the acceptable fits (at the 5% confidence level) 
there are also SB and E models, but the latter always under- 
estimate the mid-infrared fluxes. We are unable to distinguish 
between models with different formation redshifts. 

Overall, the results show that most ISOCAM cluster 
sources have SEDs typical of actively star-forming galaxies, al- 
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Table 6. Results from model SEDs for known cluster galaxies in CI 0024+1654. The columns list: LW3 identification number; 
best-fit SED model; number of data points in the observed SED; x 2 of the fit; rejection probability of the fit; other acceptable 
models (rejection probability 95% or less); 15/zm k-correction; luminosity at \5 fim (h=0.7, Q m = 0.3, Q./ = 0.7) using the 
transformation from Boselli et al. I1998K and k-corrected using best-fit SED; total infrared luminosity including the k-correction 
using the best-fit SED and computed using the transformation from Elbaz et al. J2002> ; SFR from total IR luminosity using the 
transformation of Kennicutt ( 1998): position in the V-I vs. I colour-magnitude diagram (Fig.0 where 'MS' stands for 'Main 
Sequence', 'BO' for 'Butcher-Oemler', 'U' for 'Unknown'; EW [On] line from Czoske et al. J2001> : SFR computed from the 
[O n] emission line using the transformation to the optical SFR; ratio infrared to optical SFR. 



id 


Best-fit 


(tof 


X~ 


Rejection 


Other SED 


k-corr. 




Lir 


SFR[IR] 


Colour 


EW[On] 


SFR[0 n] 


SFR[IR]/ 




model 


data points 




Probability 


models 




(L ) 


(L ) 


(Moyr- 1 ) 


Fig. 8 


(A) 


(Moyr- 1 ) 


/SFR[0 n] 


05 


S 


12 


17.8 


96.2 


none 


1.45 


2.03e+09 


9.56e+10 


16 


U 








OS 


SB+1 


15 


23.4 


97.6 


none 


1.80 


9.66e+09 


4.54e+ll 


77 


BO 








10 


S 


4 


1.9 


79.8 


S.SB+l 


1.46 


2.21e+09 


1.04e+ll 


18 


BO 








18 


s 


1 1 


3.9 


13.4 


E, S, SB+1 


1.46 


7.85e+09 


3.69e+ll 


63 


BO 








21 


s 


7 


5.4 


74.7 


SB+l.S 


1.46 


2.11e+09 


9.93e+10 


17 


BO 








22 


s 


1 1 


7.6 


52.6 


S, SB+1 


1.46 


3.96e+09 


1.87e+ll 


32 


BO 


17.5 


4.9 


6.5 


24 


s 


4 


1.2 


79.8 


E, SB+l.S 


1.42 


2.07e+09 


9.74e+10 


17 


MS 


5.0 


0.8 


21.3 


26 


s 


4 


1.4 


79.8 


S, SB+1 


1.73 


4.51e+09 


2.17e+ll 


36 










27 


SB+1 


g 


6.7 


87.8 


SB 


1.80 


5.50e+09 


2.59e+ll 


44 


MS 








Is 


S 


10 


8.6 


71.4 


E, E, S, SB+1 


1.73 


1.96e+09 


9.23e+10 


16 


BO 


9.6 


1.0 


16.0 


3s 


SB+1 


4 


1.3 


79.8 


S, SB+l.E 


1.80 


1.59e+09 


7.48e+10 


13 


MS 








4s 


E 


20 


9.4 


7.3 


E. S, SB+1 


0.78 


9.87e+08 


4.66e+10 


8 


MS 









though not necessarily observed in an exceptional starbursting 
phase. 

The SEDs of the nine cluster sources from the primary list 
were combined assuming the mean redshift of the ISOCAM 
cluster sources, and normalising the individual SEDs with the 
flux density in the rest-frame H band, either observed or as pre- 
dicted by the best-fit SED models (see also Biviano et al. 2004 1. 
We did not consider the three supplementary sources; since 
their 15 fim flux values are rather uncertain, including these 
sources would increase the scatter of the average 15^/m flux 
density value. The resulting average observed SED is shown in 
Fig- EI along with a SB+1 model fit to guide the eye (S models 
fit the average SED equally well). The average SED provides a 
synthetic view of the spectral shape of the 15 fim cluster galax- 
ies, but we refrain from a physical interpretation of it, which 
would require a detailed analysis of the systematics inherent to 
its construction. 

4.2. Infrared and optical star formation rates for cluster 
members 

The mid-infrared emission, free from dust extinction, is a re- 
liable tracer of star formation (Genzel & Cesarsky 2000). The 
infrared emission from a galaxy is the sum of various contribu- 
tions including: 

1) continuum emission from dust particles 

2) line emission from carriers of Unidentified Infrared Bands 
(UIBs) 

3) line emission from ionized interstellar gas 

4) emission from the evolved stellar population that dominates 
early-type galaxies 

5) non-thermal emission from radio sources. 

The infrared spectrum of a galaxy depends on its morpholog- 
ical type and evolutionary status. For elliptical galaxies, the 
spectrum is similar to a blackbody continuum at a temperature 



10.0 F 



z 




1 10 
rest — frame A (fim) 



Fig. 4. The average SED of 9 ISOCAM cluster sources (the 
sources from the supplementary list have not been considered). 
Error bars give the rms values for the flux density. An upper 
limit of 400 fily was used for LW2 at 7yum. The solid line 
shows a model fit (SB+1) to the average SED. 



of 4000-6000 K with a [7 jt/m]/[15 fim] flux ratio of about 4.5, 
while for spiral galaxies the ratio of the [7 //m]/[15 fim] flux is 
around 1 (Boselli et al. I1998I I. Starburst galaxies are charac- 
terized by a rapid increase in emission towards 15 fim because 
of the contribution from very small grains, which are dust par- 
ticles with radii of ~10nm that are abundant in star forming 
regions (Laurent et al. 2000 1. 

The luminosities at 15 fim, in units of solar luminosity, were 
obtained using the relationship (Boselli et al. H998i : 

Ll5^m, = 4 n D L 2 Fi 5/im 5l5^m (1) 

where F^m is the flux at 15 fim in mJy and Sis^ m = 5.04x 10 12 
Hz is the bandwidth of the LW3 filter. The total (8 - 1000 fim) 
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Fig. 5. The total infrared luminosity distribution for cluster 
members. The median value for Lfr is ~ 1.0 x 10 11 L Q 



infrared luminosity was obtained from the empirical relation of 
Elbaz et al. (120021: 



L IR = ll.l!^x(yL v [15yum]) 



0.998 



(2) 



The total IR emission was also obtained directly from the 
best-fit SED model and provides an alternative determination 
of the total infrared luminosity. The total IR luminosities ob- 
tained from the two methods were compared and found to agree 
to within 32%. We adopt the IR luminosities obtained from 
Equation [2] for the cluster galaxies. The luminosities are listed 
in Tableland include the k-corrections obtained from the best- 
fit SEDs. Depending on the best-fit model, the k-correction at 
the cluster redshift ranges from 0.8 to 1.8 with a median (mean) 
of 1.5 (1.5). The distribution of the total infrared luminosities 
for cluster galaxies is given in Fig.|5] The mid-infrared cluster 
members have total IR luminosities between 4.7 x 10 10 L and 
4.5 x 10" L , with a median value of 1.0 x 10 11 L . Six of 
the 12 galaxies have total infrared luminosities above 10 11 L , 
which classify them as LIRGs (10 11 L < Lir < 10 12 L , e.g. 
Genzel & Cesarsky 2000 1, and four more have infrared lumi- 
nosities above 9 x 10 10 L . The star formation rates in units of 
solar masses per year, SFR[IR], were derived using the relation 
of (Kennicutt 1998 1: 



SFR[IR] 



1.71xl0- ll) (L IR /L o ) 



(3) 



The values of SFR[IR] are listed in Table[6] The infrared SFRs 
range between 8 M Q yr _I and =2 77 M© yr _I , with a median 
(mean) value of 18 (30) M yr 1 . 

Half of the 15/im cluster galaxies are LIRGs and four 
more are within 1 <x of the luminosity of a LIRG. Star form- 
ing episodes that are enshrouded by dust are commonly found 



in these galaxies in the field. The cause of the starburst phase 
of LIRGs in the local universe is still debated, but in the range 
of infrared luminosities measured with ISOCAM the starburst 
seems to be linked to mergers or interactions between pairs of 
spiral galaxies of unequal size and single galaxies in about one- 
quarter of the sources (Hwang et al. 1999; Ishida & Sanders 
2001 1. The cause of the starburst in the 15//m sources of CI 
0024+1654 seems to be consistent with what is observed in 
the local universe because of the evidence for interactions and 
mergers in the HST maps (see Sect. 15. 2> . 

The [0 11] line is often used to determine the SFR in galax- 
ies at intermediate redshift, when the Balmer hydrogen lines 
are outside the spectral range (e.g. Barbara & Poggianti ll9971 
Jansen et al. 2000 1. There are several problems associated with 
this line. [On] is subject to dust extinction, and its strength de- 
pends heavily on the metallicity and ionization of the interstel- 
lar medium (Kennicutt 1998 Jansen et al. 2000). Furthermore, 
Balmer lines are more correlated to UV emission from young 
stars than [On]. Therefore, the uncertainties in SFR computed 
from [On] are very large. Nevertheless, there may be times 
when it is the only indicator available, in which case finding 
ways to correct the determination for dust extinction will be 
important. It may well yield only a lower limit on the SFR. 

The Equivalent Width (EW) of the [O n] line was available 
from Czoske et al. (2001 ) for only three of the mid-infrared 
cluster members. The [On] line is redshifted into the V band 
and can be used to compute the optical SFR provided that the 
V-band magnitude is known. The luminosity L v was derived 
from the V-band magnitude in units of erg s _I A -1 using: 



L v = 4 n D 2 x (3.08 x 10 24 ) 2 x l(T a4m " x 3.92 x 10~ 9 



(4) 



where Dl is in Mpc and m v is the apparent magnitude in the 
V-band. The luminosity of the [On] line is given by L[On] 
= EW[On]x Ly in units of erg s _1 . Finally, the SFR was ob- 
tained using SFR[0 n] ~ 1 .4 x 10~ 41 L[0 n], in units of M yr _ 1 
(Kennicutt H9^8l . The SFRs derived from the EWs of the [0 11] 
lines are listed in Table [6] and should be regarded as lower 
limits because no correction for extinction was applied. The 
mean value of the optical SFR is 2.2 M yr~' and increases to 
3 M yr~' when the canonical value of 1 mag at Ha (Kennicutt 
1992 1, for extinction in the optical, is applied. The values of 
SFR[On] show a wide range, and are sometimes more than 1 
order of magnitude lower than those obtained in the infrared. 
Recent data show that there is a better agreement when the 
SFRs are computed from the Ho- line and the mid-infrared 
(Kodama et al. 120051 The ratios of SFR[IR]/SFR[0 n] range 
from 7 to 21 with a median value of 16 (section l4~2l . Therefore 
the vast bulk of the star formation is missed when the [O n] line 
emission is used. A large fraction of the star formation in CI 
0024+1654 is enshrouded by dust (see Cardiel et al. l2003l for 
a thorough discussion on the comparison of SFR-estimators). 
Due et al. ( 2002 1 and Balogh et al ( 2002 1 have arrived at a very 
similar conclusion for the cluster Abell 1689, which has many 
mid-infrared sources. Similar results have also been obtained 
for field galaxy samples (see, e.g., Chariot et al. 2002 1. 

In our sample, AGN activity does not seem to be a major 
contributor to the energy output of the systems. Measurements 
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Fig. 6. (a) Adaptive Kernel map for cluster members. Black boxes represent the 15 fim sources while dots and isodensity contours 
trace the distribution of component-A galaxies, (b) Radial distribution from the centre of component-A galaxies (continuous line) 
and 15//m cluster galaxies (dashed line) within the region covered by the ISOCAM observations. 



of the Ho- and [N n] emission lines are required to further de- 
termine the AGN contribution to the mid-infrared emission 
(Balogh et al. 2002 1 and future radio observations could also 
help in determining the AGN contribution (Rengarajan et al. 
rrW7l Dwarakanath & Owen fT999l . 

5. Spatial, dynamical and colour properties of the 
15//m cluster galaxies 

5.1. Spatial and velocity distributions 

The spatial and radial distributions of the cluster galaxies, with 
and without detectable 15 fim emission, are given in Fig.il The 
isodensity contours and dots in Fig. |6(a)| refer to galaxies be- 
longing to component A, while the black boxes represent the 
cluster members with 15/im emission. The figure shows that 
the cluster members detected at 15 /mi are less concentrated 
than the cluster galaxies not detected in the mid-infrared. The 
Kolmogorov-Smirnov test reveals that there is only a 4% prob- 
ability that the two distributions are drawn from the same par- 
ent one. 

The velocity distributions for cluster galaxies with and 
without detectable 15 fim emission are given in Fig. |2] The 
counterparts of the 15 fim sources have redshifts which place 
them in Component A - the larger of the two interacting com- 
ponents. One source, ISO_C10024_12, with z = 0.386 is close 
to the boundary. The velocity dispersion computed in the clus- 
ter rest-frame (Fig fj}, using the biweight estimator (Beers et 
al.[T990l in the central region of the ISO map, is cr v = 914+^ 



km s for a sample of galaxies, while it is <x v = 979*™ km 
s _1 for the 15/im galaxies. The 15 fim sources are associated 
mainly with spiral and emission-line galaxies, and studies of 
nearby clusters indicate a larger velocity distribution for this 
population relative to the general cluster population (Biviano 
et al. nWfl . 

5.2. Colour-magnitude diagramme 

The colour-magnitude diagramme for cluster galaxies is given 
in Fig. [7] It provides an important link between the optical 
and mid-infrared properties of the ISOCAM galaxies. Apparent 
magnitudes in the V and I bands are available from Czoske et 
al. i200Tl for 11 of the 13 ISOCAM cluster galaxies. The cir- 
cles represent component-A galaxies within the region mapped 
by ISOCAM. The spectroscopically confirmed mid-infrared 
cluster galaxies are indicated by filled circles and are labelled 
with the ISO source number. Section M of Fig.Qcontains the 
galaxies that satisfy the BO definition, given in section ^ Six 
of the 15/im galaxies fall in section M and satisfy the BO re- 
quirement. Only 20% of the cluster galaxies in section M are 
detected at 15/im. However, this fraction increases to ~38% 
when only galaxies brighter than mi = 19.8 are considered. 
The non-detections of fainter I-band sources are probably due 
to the sources falling below the sensitivity limits of ISOCAM. 

All 15/mi cluster sources have in common an excess of 
mid-infrared emission that sets them apart from the rest of the 
cluster members. In at least some of these systems, one would 
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Fig. 7. V - I colour-magnitude diagramme for component-A galaxies (empty circles) within the region mapped by ISOCAM. 
The eleven 15 fim cluster galaxies (filled circles) are identified by their ISO source number from Tables |2] and [3] The V and I 
magnitudes were not available for source ISO_C10024_12. The galaxies in section M, which is defined by V - 1 < 1.8 and I < 
20.4, are BO galaxies. Surprisingly, five of the 15 fim sources are in section N and not associated with BO galaxies. 



expect that the star forming activity is primarily triggered by 
interactions and mergers. Recent changes in the galaxy popula- 
tion reveal themselves in the mid-infrared and help identify the 
processes that cause the burst of obscured star formation. The 
infall of smaller groups into the cluster environment provides 
a way of promoting slow encounters and mergers within clus- 
ters (Mihos 2003 1. Slow encounters between galaxies are better 
able to drive instabilities than fast encounters. Complex interac- 
tions of the kind that may occur between galaxies in this cluster 
have been studied in the local Universe in Hickson Compact 
Groups (see, e.g., Lopez-Sanchez et al. 2004 and references 
therein). 

We looked for indications of ongoing and past interaction 
events among the 15 /mi cluster sources, either in their mor- 
phologies, or in their spectra. A HST image was available for 
part of CI 0024+1654 from the MORPHS collaboration (Smail 
et al. 11997b . This image is limited to the South-Eastern re- 
gion of the cluster core and contains seven of the 15/im clus- 
ter sources. For these galaxies, miniature maps centred on the 
ISOCAM source coordinates and having an area of 30 x 30 
arcsec 2 were obtained. The miniature maps are presented in 
Fig. El 

Miniature maps are available for 3 of the BO galaxies in 
section M of Fig. 13 The HST maps show that ISO_C10024_18 
and ISO_C10024_22 do not have companions within a radius of 
5" while ISO_C10024_08 appears to be interacting with at least 
another galaxy. 

HST miniature maps for four of the non-BO 15/im emit- 
ters in section N are also given in Fig. [8] and each has at least 
one nearby companion. ISO_C10024_4s appears to be merging 
because it has multiple nuclei and has at least two other com- 
panions. ISO_C10024_24 and ISO_C10024_3s show evidence of 



ongoing interactions (e.g. tidal bridges) with at least one nearby 
galaxy. 

Even though the number of infrared sources with HST 
maps is small, there seems to be a difference between the BO 
and non-BO galaxies in the sense that the former tend not 
to have nearby companions whereas the latter do. The op- 
tical counterparts of the 15/im sources on the main colour- 
magnitude sequence are not passive early type galaxies as one 
would expect from their position in the colour-magnitude dia- 
gramme. Rather, the optical counterparts appear to be galaxies 
that are involved in interactions and mergers in the HST im- 
ages. 

There is indirect support for this apparent difference be- 
tween the BO and non-BO galaxies from the results of the SED 
fitting process. Five of the six 15 fim galaxies in the BO region 
of the diagram have SEDs that are best-fit by S-type models. 
On the other hand, two out of four 15 fim galaxies on the main 
sequence of the colour-magnitude diagram have SEDs that are 
best-fit by SB +1 -type. If the starburst events are triggered by 
interactions, there is marginal evidence from the SED analysis 
that the 15 fim galaxies on the main colour-magnitude sequence 
have suffered more interaction events in the recent past than the 
BO 15 fim galaxies. 

Given the scarce morphological data, and the fact that the 
SED fitting results are not very well constrained (SB+l-type 
models also provide acceptable fits to the SEDs of BO galaxies, 
see Table |6j, we consider the evidence for a difference in the 
interaction properties of BO and non-BO 15/im galaxies only 
tentative. We conclude from our data that interaction events do 
trigger the star formation activity (and stimulate the IR emis- 
sion) in at least some 15 fim cluster galaxies. 
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Fig. 8. Miniature maps obtained from HST. The maps are centred on the 15 fim source coordinates and labelled with the source 
number in Tables [2] and [3] The letters M and N at the top right of the images give the positions of the cluster sources in the 
colour-magnitude diagramme (Fig. [7}. All images are 30 x 30 arcsec 2 . Insets show details of the optical counterpart. North is up, 
East to the left. 



Obscuration by dust could explain why some 15 /mi sources 
lie on the cluster main sequence. We tested this possibility by 
computing the de-reddened V and I magnitudes of the 15/im 
galaxies, using the dust-free best-fit model SEDs, as given by 
the GRASIL code. As expected, once de-reddened the 15/im 
galaxies on the colour-magnitude sequence move in the colour- 
magnitude diagramme toward the BO region, and in one case 
(ISO_C10024_3s), it even becomes part of the BO category 
(Fig-0- Dust obscuration at least partially explains the unusual 
colours of the 15 /mi sources on the cluster main sequence. 

6. Comparison of CI 0024+1654 with other clusters 
observed by ISO 

Observations of clusters of galaxies at high redshift (z > 0.4) 
with SCUBA have yielded an excess of sub-millimeter sources 
with high values of the SFR (Best 2002). Here we compare the 
results of the observations of CI 0024+1654 with four clusters 
of galaxies that were observed with ISOCAM at 15 //m (Fadda 
et al. l2000l Metcalfe et al. l2003l to determine if C10024+1654 
has more LIRGs than the other clusters. We have not included 
the cluster J 1888.16 CI in Table0because the relevant data 
have not been published (Due et al. 2004). However this clus- 
ter seems to be comparable to C10024+1654 because it has 
many infrared sources and at least 6 of them are confirmed 
LIRGs. The characteristics of the observations, and results for 
the five clusters, are summarized in Table The observed ar- 
eas and sensitivities are comparable for the first three clusters 
in Tabled but the area is somewhat smaller for Abell 2218, and 
much smaller for the ultra-deep observations of Abell 2390. 



The number of 15 /mi sources that are identified with clus- 
ter galaxies ranges from 13 in CI 0024+1654 to 1 in Abell 370. 
The number of 15 /mi sources that have fluxes consistent with 
LIRGs, within the precision of the measurements is listed in 
Table □ In CI 0024+1654 the number of LIRGs is 10 (Table 
0) and includes 6 sources with L/« > 10 11 L Q and 4 more with 
L/r between 9 x 10 10 L G and 1 x 10 11 L G (Table|6}. The number 
of LIRGs detected in the other clusters is either or 1 (Table 
13). The virial mass and radius of CI 0024+1654, Abell 370 and 
Abell 1689 are the same to within 20% and are considerably 
smaller than the values for the more massive clusters Abell 
2218 and Abell 2390. 

However the ratio of the virial mass to area within the virial 
radius varies by less than 50% for the 5 clusters. We now com- 
pare the number of LIRGs in CI 0024+1654 with those in the 
other clusters. The comparison is limited to LIRGs because 
it eliminates selection effects for the 15/im sources. In any 
case most of the 15 /mi sources in the nearer clusters are too 
faint to be detected in the more distant clusters CI 0024+1654 
and Abell 370. In the comparison, the number of LIRGs in CI 
0024+1654 was multiplied by the ratios of : 

a) the virial mass per unit area (since cluster mass is approxi- 
mately proportioned to richness - see Bahcall & Cen 1993 ) 
of the cluster to that of CI 0024+1654 

b) the square of the distance to the cluster and that to CI 
0024+1654 

c) the observed solid angle of the cluster to that of CI 
0024+1654. 
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Table 7. Summary of ISOCAM observations and results at 15 /mi for five clusters of galaxies. The data was obtained from 
Metcalfe et al. J20l)3l for Abell 370, Abell 2218 and Abell 2390, Fadda et al. ( 120001 and Due et al. ( 120021 for Abell 1689, and 
this paper for CI 0024+1654. The content of the columns in the table are as follows: name and redshift of the cluster, total area 
scanned, sensitivity reported at the 5<x level, flux of the weakest reported source in //Jy, total observation time, total number of 
sources detected including sources without redshifts and stars. Then number of cluster galaxies, virial radius, virial mass, number 
of sources with Lir > 9 x 10 10 L0 detected and expected. The expected number of sources was obtained by comparison with CI 
0024+1654 as described in the text. Virial radii and masses are from Girardi & Mezzetti ( 1200 1 1 > and King et al. d2002t . 



Cluster 


z 


Area 


Sensitivity 


Faintest source 


Obs. t. 


Tot. n. 


Cluster 


Rvir 


M vlr 


n. of IR 


sources 






C 2 ) 


(5o-) 


G"Jy) 


(sec) 


sources 


galaxies 


(lr 1 Mpc) 


(lr 1 10 14 M ) 


Detected 


Expected 


C10024 


0.39 


37.8 


140 


141 


22615 


35 


13 


0.94 


6.42 


10 




A370 


0.37 


40.5 


350 


208 


22688 


20 


1 


0.91 


5.53 


1 


8 


A1689 


0.18 


36.0 


450 


320 


9500 


18 


11 


1.1 


5.7 





1 


A2390 


0.23 


7.0 


100 


54 


29300 


28 


4 


1.62 


20.35 





1 


A2218 


0.18 


20.5 


125 


90 


22000 


46 


6 


1.63 


18.27 





1 



The values are listed in the last column of Table0 The simple 
scaling of the mapped area will not be sufficient if the infrared 
galaxies have a spatial distribution that is different from the 
overall cluster population. In particular the mapped regions of 
the three nearer clusters (AbeU 1689, Abell 2218, Abell 2390) 
are smaller and more confined to the central regions than the 
two more distant clusters (Abell 370 and CI 0024+1654). The 
infrared sources could preferentially occur in the peripheral re- 
gions of clusters where dynamical interactions between galax- 
ies may be more effective in triggering interactions and bursts 
of star formation (Mihos 2003). In this case, the LIRG pop- 
ulations for the complete clusters may in some cases be less 
different than the regions currently mapped would make it ap- 
pear. Future observations with Spitzer are needed to address 
this possibility (Werner et al. 2004 1. 



6.1. LIRGs in CI 0024+1654 and Abell 370 



The clusters CI 0024+1654 and Abell 370 have comparable 
redshifts and were observed by ISOCAM with almost identical 
observational parameters, such as, on-chip integration time, to- 
tal area observed in raster mode, and total observation time. 
The total number of cluster sources in CI 0024+1654 is 13 
whereas the corresponding number is 1 for Abell 370 (Table0. 
The observed number of LIRGs in Abell 370 is only 1 whereas 
8 were expected from the comparison with CI 0024+1654. 
There is a real difference between the populations of LIRGs 
and also in the ratio of the mass to infrared light for the two 
clusters. No LIRGs were detected in the three clusters Abell 
1689, AbeU 2218 and Abell 2390 at lower redshifts whereas a 
total number of 3 was expected from the comparison with CI 
0024+1654. A population of luminous infrared galaxies is also 
absent from the central regions of the Coma and Virgo clus- 
ters (Boselli et al. fT997lfT998l Q uillen et al. [T999l Leech et al. 
1999| Tuffs et al.Ep02a 2002b|i. 



6.2. Mid-infrared sources in Abell 1689 and Abell 2218 

It is, further, interesting to compare the two clusters Abell 1689 
and Abell 2218 that are at the same redshift. The observations 
of Abell 2218 are more sensitive to 15 /mi sources than those of 
Abell 1689 (Table0. However there is no significant difference 
between the total number of 15 /mi sources when allowance is 
made for the scanned area and virial mass per unit area. The dif- 
ference caused by the sensitivity of the observations is revealed 
by a comparison of the median fluxes of the 15//m sources 
which is ~ 600 juJy for Abell 1689 and ~ 150 /Jy for Abell 
2218. All of the 15 /mi sources in Abell 1689 would have been 
easily detected in Abell 2218, if present, whereas only one of 
the sources in Abell 2218 might have been detected in the Abell 
1689 measurement. The 15 fim sources detected in Abell 2218 
have lower fluxes and luminosities than those in Abell 1689 and 
this effect is not caused by the sensitivities of the observations 
of the two clusters. It is very interesting that Abell 1689 and 
Abell 2218 follow the same trend identified in the comparison 
between CI 0024+ 1 654 and Abell 370. 

Large numbers of cluster sources are detected by ISOCAM 
at 7 /mi in both Abell 1689 and Abell 2218. In Abell 2218 the 
SEDs of most of these 7 /mi sources are well fit by models of 
quiescent ellipticals with negligible SFRs and a median lumi- 
nosity of ~ 8 x 10 8 L (Biviano et al. l2004t . 

6.3. Dynamical status of clusters and population of 
mid-infrared sources 

CI 0024+1654 is involved in a recent merger, as shown by the 
two velocity components in Fig. [5] Abell 1689 also has a com- 
plex and broad velocity structure, with evidence for three dis- 
tinct groups that overlap spatially and are well separated in 
velocity space (Girardi et al. I1997I I. The two clusters Abell 
370 and Abell 2218 show no evidence for major merger ac- 
tivity. The number of luminous 15 /mi sources can be, in part, 
explained by considering the dynamical status of the clus- 
ter. Simulations made by Bekki (1999 1 show that the time- 
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dependent tidal gravitational field existing in cluster-group 
mergers induces secondary starbursts by efficiently transferring 
large amounts of gas from the disk to the nucleus. The model 
establishes a link between the population of starburst or post- 
starburst galaxies and the presence of substructures in clusters. 
The existence of starburst and post-starburst galaxies in a merg- 
ing cluster, and spread over a wide area, is an important predic- 
tion of this model which seems to be in agreement with the 
results from CI 0024+1654 and Abell 1689. 

Furthermore it can be expected that merging clusters at 
much larger redshifts than CI 0024+1654 will contain LIRGs 
and ULIRGs because the interacting and merging galaxies will 
be more gas rich (see, e.g. Balogh et al. 2000 Kauffmann & 
Haehnelt l2000l Francis et al. l200ll . 

In this context it is interesting to note that Dwarakanath & 
Owen ( 1999 1 found different radio source populations in two 
very similar clusters. The number density of radio sources in 
Abell 2125 (z = 0.246) exceeded that in Abell 2645 (z = 0.25) 
by almost an order of magnitude. The cluster Abell 2125, with 
the larger number of radio sources, also shows evidence for a 
merger. 

6.4. Mid-infrared sources in CI 0024+1654 and Abell 
1689 

Many infrared sources were found in the cluster Abell 1689 
(Fadda et al. 2000 1 that is at a much smaller redshift than CI 
0024+1654. No galaxy in Abell 1689 has total infrared lumi- 
nosity above 1 x 10 1 1 L (we use our cosmology and Equation|2] 
to convert the 15 /mi fluxes from Fadda et al. (2000) to total 
infrared luminosities), while half of the 15 /mi cluster sources 
in CI 0024+1654 are above this luminosity. All of the 15 /mi 
sources detected in CI 0024+1654 would have been easily de- 
tected in Abell 1689, if present, whereas most of the sources 
in Abell 1689 would not have been detected in the observa- 
tions of CI 0024+1654 when allowance is made for different 
sensitivities and distances. The two clusters have very similar 
virial radii and masses (Table 0. However a much larger part 
of the periphery of CI 0024+1654 was observed because it is 
at a higher redshift. It is possible that some of the differences 
between the 15 //m sources arises from a change in the source 
population between the core and periphery of the two clusters. 

The population of 15 /mi sources may also be influenced by 
the history of the cluster. Various processes such as ram pres- 
sure stripping of gas from cluster galaxies, tidal effects and pre- 
vious starbursts will inevitably leave less gas in the interacting 
and merging galaxies to fuel the starburst (e.g. Gunn & Gott 
\W7H Byrd & Valtonen [T33U| Fuiita fT9^8l l. These processes 
may vary from one cluster to another depending on its history 
and provide dispersion in the luminosity of the starbursts. The 
large number of LIRGs in CI 0024+1654 are fuelled by the gas 
rich progenitor galaxies. It is conceivable that in Abell 1689 the 
starbursts are less luminous either because there is a smaller 
supply of gas to fuel the outburst, or the luminosity has de- 
creased because a longer time has elapsed since the last merger 
or there is a significant difference between the sources in the 
core and periphery. 



The average value of the ratio SFR[IR]/SFR[0 n] is rea- 
sonably comparable for both clusters (it is ~ 15 for CI 
0024+1654 and ~ 11 for Abell 1689). The internal proper- 
ties of the 15 /mi sources in CI 0024+1654 and Abell 1689 are 
broadly similar even though their distributions of the luminosi- 
ties of 15/mi sources are different. The colour-magnitude dia- 
grammes of the 15 //m sources in the two clusters are also quite 
similar. 

7. Conclusions 

The cluster CI 0024+1654 was observed with ISO. A total of 
35 sources were detected at 15/mi and all have optical coun- 
terparts. Sources with known redshift include four stars, one 
quasar, three background galaxies, one foreground galaxy and 
thirteen cluster galaxies. The remaining 13 sources are likely 
to be background sources lensed by the cluster. 

The spatial, radial and velocity distributions were obtained 
for the cluster galaxies. The ISOCAM cluster galaxies appear 
to be less centrally grouped (in the cluster) than those not de- 
tected at 15 fim and the Kolmogorov-Smirnov test reveals that 
there is only a 4% probability that the two distributions are 
drawn from the same parent population. No statistically signif- 
icant differences were found between the velocity distributions 
of the 15//m sources and other cluster galaxies in the region 
mapped by ISOCAM. 

Spectral energy distributions were obtained for cluster 
members and used as indicators of both morphological type 
and star forming activity. The ISOCAM sources have as best- 
fit SEDs predominantly those of spiral or starburst models ob- 
served 1 Gyr after the main starburst. Star formation rates were 
computed from the infrared and the optical data. The SFRs in- 
ferred from the infrared are one to two orders of magnitude 
higher than those computed from the [O n] line emission, sug- 
gesting that most of the star forming activity is hidden by dust. 

A colour-magnitude diagramme is given for cluster sources 
falling within the region mapped by ISOCAM. V and I-band 
magnitudes are available for 1 1 of the cluster sources, and 6 of 
these have colour properties that are consistent with Butcher- 
Oemler galaxies and best-fit SEDs that are typical of spiral 
models. The remaining 15/mi cluster galaxies have colours 
that are not compatible with Butcher-Oemler galaxies and have 
best-fit SEDs that are typical of starburst galaxies 1 Gyr after 
the main burst. HST images are available for these latter sys- 
tems and all have nearby companion galaxies. These results 
suggest that interactions and mergers are responsible for some 
of the luminous infrared sources in the cluster. 

The 15/mi sources in CI 0024+1654 were compared with 
four other clusters observed with ISOCAM. The results show 
that the number of LIRGs in Abell 370 is smaller than expected 
by about one order of magnitude, if Abell 370 were to be com- 
parable with CI 0024+1654. Furthermore no LIRGs were de- 
tected in Abell 1689, Abell 2218 and Abell 2390 when a total 
of 3 was expected, based on the results from CI 0024+1654. 
A comparison of the mid-infrared sources in Abell 1689 and 
Abell 2218 shows that the sources in Abell 1689 are more lu- 
minous than in Abell 2218 and follow the same trend identified 
in the comparison between CI 0024+1654 and Abell 370. There 
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is clear evidence for an ongoing merger in CI 0024+1654 and 
Abell 1689. The number and luminosities of the mid-infrared 
cluster sources seem to be related to the dynamical status and 
history of the clusters. 
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Appendix A: Additional notes on some 
mid-infrared sources 

ISO_C10 024_02 Radio quiet quasar PC 0023+1653 (Schmidt 
et al. 119861 at a redshift of 0.959 with an X-ray luminosity 
L x = 1.4 x 10 44 ergs-' (Soucail et al. l2000l 

ISO_C10024_04 The 15 /mi map has two bright regions that 
may contain contributions from two optical counterparts. 
There are no measured redshifts. 

ISO_C10024_06 The 15/mi emission appears to have contri- 
butions from various galaxies in the field-of-view. No red- 
shifts are available for these sources. 

ISO_C10024_10 Most of the 15 /mi emission is centred on a 
cluster galaxy at redshift z = 0.400 (Table the bright- 
est of the two galaxies in Fig. [8] The model fit to the SED 
is an Sc (Fig. The source was also detected in X-rays 
with ROSAT (Bohringer et al. 120001 The optical identifi- 
cation for the ROSAT source is uncertain because there are 
two sources in the error box. One of the sources is a typi- 
cal cluster galaxy while the second is a star-forming fore- 
ground galaxy (z = 0.2132). The cluster galaxy is closer to 
the 15 /mi coordinates and was adopted as the mid-infrared 
counterpart. It is the only cluster AGN detected so far in CI 
0024+1654. 

ISO_C10024_16 The redshift of the optical counterpart on the 
VLT image is not known. The 15 fim emission is centred on 
one galaxy and partially extends to two very faint optical 
sources that are north of the main optical source. 

ISO_C10024_17 The redshift of the optical counterpart on the 
VLT image has not been measured. The 15 fim emission is 
centred on the bright optical source and is partially elon- 
gated towards two other sources on the VLT map. 

ISO_C10024_18 The optical counterpart is a large, edge-on, 
late type, cluster spiral galaxy. The 15 /mi emission is cen- 
tred on the galaxy and has contributions from several other 
sources. The HST image reveals a disturbance of the spiral 
structure. 

ISO_C10024J22 The optical counterpart is a cluster galaxy 
with z = 0.3935. The HST image (Fig. |SJ reveals a spi- 
ral galaxy with an inner ring, a smoother outer arm (Smail 
et al l 19971 1 and bright knots. 

ISO_C10024J23 The redshift of the optical counterpart is not 
available. At the 15 //m coordinates, the HST image shows 
two interacting spiral galaxies and a tidal arm. 



ISO_C10024_30 The source is outside the boundaries of the 
VLT and HST images (Fig.0 an d has a faint optical coun- 
terpart (Fig. 4d in Czoske et al. l2001t . 
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